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We perform a lattice study of charmonium-like mesons with J pc — 1 ++ and three quark contents 
ccdu, cc(uu + dd) and ccss, where the later two can mix with cc. This simulation with Nf = 2 and 
m-n ~ 266 MeV aims at the possible signatures of four-quark exotic states. We utilize a large basis 
of cc, two-meson and diquark-antidiquark interpolating fields, with diquarks in both antitriplet and 
sextet color representations. A lattice candidate for A' (3872) with I = 0 is observed very close to the 
experimental state only if both cc and DD * interpolators are included; the candidate is not found 
if diquark-antidiquark and 7272* are used in the absence of cc. No candidate for neutral or charged 
A(3872), or any other exotic candidates are found in the 7=1 channel. We also do not find signatures 
of exotic ccss candidates below 4.2 GeV, such as Y (4140). Possible physics and methodology related 
reasons for that are discussed. Along the way, we present the diquark-antidiquark operators as linear 
combinations of the two-meson operators via the Fierz transformations. 


I. INTRODUCTION 

The experimental discovery of charged resonances 
27 c (3900) + [I[ and Z(443[2, 0 gives signatures for 
hadrons with minimal quark content ccdu. The neutral 
A(3872) and yet-unconfirmed Y(4140) with charge par¬ 
ity C = +1 also appear to have significant four-quark 
Fock components. Most of the observed exotic states 
have J p = 1 + . The J p for some has not been settled ex¬ 
perimentally and J p = 1 + presents one possible option. 

In this paper, we perform a lattice investigation of 
the charmonium spectrum, looking for charmonium-like 
states with quantum numbers J PC = l +_,_ and three 
quark contents: ccdu , cc(uu + dd) and ccss , where the 
later two channels have 7 = 0 and can mix with cc (C 
indicates C-parity of neutral isospin partners for charged 
states). Our main interest in these channels is aimed at 
a first-principle study of A (3872) and Y(4140), which 
were observed in A(3872) — > J/ifp, J/ifoj, DD* and 
Y(4140) -A J/if<t>, for example. 

From the experimental side, the long known exotic 
candidate A(3872) Q is confirmed to have J PC = 1 ++ 
[S:]. However, questions about its isospin remain unset¬ 
tled. If it has isospin 7=1, one expects charged part¬ 
ners. Observation of a nearly equal branching fraction 
for 27(3872) —> J/ifu and 27(3872) -4 J/ifp decays @ 
and searches for charged partner 27(3872) states decay¬ 
ing to [Z] speak against a pure 7 = 1 state. There 

are a few other candidates with C = +l that couldpossi- 
bly have J PC = 1 ++ like 27(3940) Q, Y(4050) ± [sj and 
27(4250) ± 0. A detailed review on these can be found 
in Ref. [§|. 
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The growing evidence for the Y (4140) resonance in the 
J /^0 invariant mass 0 serves as promising signature 
for exotic hadrons with hidden strangeness. Similarities 
in the properties of 27(3930) and Y(4140) led to an in¬ 
terpretation that 27(3930) may be a 72*72* molecule and 
y(4140) is its hidden strange counterpart 72*72* molecule 
11]. However, the upper limit for the production of 
Y(4140) in 77 —> J/ipcf is observed to be much lower 
than theoretical expectations for a 72*72* molecule with 
J PC = 0++ and 2 ++ 0]. Hence the quantum num¬ 
bers of Y (4140) stay unsettled and it remains open for a 
J pc = 1 ++ assignment. 

From a theoretical perspective, the description of such 
resonances is not settled. Several suggestions have been 
made interpreting them as mesonic molecules fl3| , as 
diquark-antidiquark structures 14], as a cusp phenomena 
[la ] or as a | ccg) hybrid meson [16j. A great deal of theo¬ 
retical studies are based on phenomenological approaches 
like quark model, (unitarized) effective field theory and 
QCD sum rules (see reviews Q). 

It is paramount to establish whether QCD supports 
the existence of resonances with exotic character using 
first principles techniques such as lattice QCD. Simu¬ 
lations that considered only cc interpolators could not 
provide evidence for A(3872). The first evidence from a 
lattice simulation for 27(3872) with 7 = 0 was reported 
in Ref. 0 , where a combination of cc as well as 7272* 
and J/ifui interpolators was used. Recently, another cal¬ 
culation using the Highly Improved Staggered Quark ac¬ 
tion also gave evidence for A (3872), using cc and 7272* 
interpolating fields 0]. The search for the Y(4140) reso¬ 
nance was performed only in [l£j ] , where a phase shift for 
J/if 4> scattering in s-wave and p -wave was extracted from 
Nf = 2 + 1 simulation using twisted boundary conditions, 
and neglecting strange-quark annihilation. The resulting 
phase shifts did not support existence of a resonance. 

The novel feature of the present study is to add 
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diquark-antidiquark [cq]g[cq)g operators to the basis of 
interpolating fields and to extend the extraction of the 
charmonium spectrum with J PC = 1 ++ to a higher 
energy range. This is the first dynamical lattice cal¬ 
culation involving diquark-antidiquark operators along 
with several two-meson and cc kind of interpolators to 
study A(3872) and Y(4140). We consider the color 
structures Q = 3 C ,6 C for diquarks, which have been 
suggested already in the late seventies [20]. Recently 
many phenomenological studies Bill and a few lat¬ 
tice studies ii| used them to extract the light and 
heavy meson spectra. In Ref. [23| a calculation us¬ 
ing two-meson and diquark-antidiquark interpolators was 
performed to investigate mass spectrum of 1 ++ exotic 
mesons in quenched lattice QCD. However, only one en¬ 
ergy level was extracted, which is not sufficient to provide 
evidence for A(3872) or Y(4140). 

In this paper we address the following questions: Is 
the lattice candidate for 77(3872) reproduced in presence 
of diquark-antidiquark operators? Which are the crucial 
operator structures for its emergence? How important 
are the [cg]g[c< 7 ]g Fock components in the established 
77(3872)? Do we find a lattice candidate for charged 
or neutral 77(3872) with 7 = 1? Do operators with hid¬ 
den strangeness render a candidate for Y(4140)? Do we 
find candidate states for other possible exotic states in 
the channels being probed? 

The paper is organized as follows. Section ITTl addresses 
the expected two-meson scattering channels below 4.2 
GeV. The lattice methodology is discussed in Sect. IIIII 
In Sect. m and the Appendix we discuss the relations 
between our diquark-antidiquark and two-meson interpo¬ 
lators via Fierz transformations. Section El is dedicated 
to results and we conclude in Sect. eh 

II. TWO PARTICLE STATES IN LATTICE QCD 

A major hurdle in excited-state spectroscopy is that 
most of the states lie above various thresholds and decay 
strongly in experiments. All states carrying the same 
quantum numbers, including the single-particle and mul- 
tiparticle states, in principle contribute to the eigenstates 
of the Hamiltonian. The determination of scattering 
properties relies on precise identification of all the eigen¬ 
states below and close above the energy of our interest. 
The continuous spectrum of scattering states in the con¬ 
tinuum gets reduced to a discrete set of eigenstates, be¬ 
cause lattice momenta are discretized due to the finite 
lattice size. 

Considering two-meson states with total momentum 
zero and without interaction, their energies are just the 
sum of the individual particle energies 

-®Ml(n)M 2 (-n) = E l(p)+ E 2 (jp) , p = — , n G N 3 . (1) 

In the presence of interactions, the energies get shifted 
depending on the interaction strength. For our lattice 


setting the noninteracting two-meson levels with J PC = 
1 ++ , total momentum zero in the indicated energy range 
are 


• 1 = 0; cc{uu + dd) and cc; E < 4.2 GeV 

D(0)D*(0), JM0M0), D(1)D*(-1), 

J/^( iM-i), rj c (l)cj(- 1), Xci(0)cr(0) . 

• 7 = 1; ccdu ; E < 4.2 GeV 

D(0)D*(0), JMO)p(O), D(1)D*(—1), 

^/V ; (l)p(-1), Xci(lM—1), Xco(l)7r(—1) . 

• 7 = 0; ccss and cc; E < 4.3 GeV 

d s (o)d* s (o), j/mm , D a (i)D* a (-i), 

JM!)</>(-!) ■ 


The parentheses denote meson momenta in units of 27 t/L. 

We consider the flavor sectors cc{uu + dd) and ccss 
separately. In nature these two 7 = 0 sectors can mix 
and they could in principle mix also in our simulation 
without dynamical strange quarks. However, if both fla¬ 
vor sectors would be treated together, then 6 + 4 = 10 
two-particle 7 = 0 states are expected below 4.2 GeV. 
This would make the resulting spectrum denser and nois¬ 
ier, so the identification of eigenstates and the search for 
exotics would be even more challenging. We therefore 
consider these two sectors separately in this first search 
for possible exotics in the extended energy region. The 
corresponding assumptions will be discussed for each fla¬ 
vor channel along with the results. 

The noninteracting energies will be shown by the hor¬ 
izontal lines in our plots, and follow from the masses 
and the single meson energies determined on the same 
set of gauge configurations [24l-(26| . The energies of the 
cr meson using single-hadron approximation are am a = 
0.302(15) and a 7? cr ( 1 ) = 0.534(22). Including two-meson 
operators up to 4.2 GeV at m n = 266 MeV should be suf¬ 
ficient in searching for narrow exotic candidates below 4.2 
GeV. Details of all the interpolators used, including the 
diquark-antidiquark interpolators, can be found in the 
next section. 

The mesons R = p, a are resonances that decay to 
7 T7r or irr] in QCD with Nf = 2. A proper simulation 
which would consider the three-meson system f2?J has 
not been performed in practice yet. In absence of this, a 
simplifying approximation for channels containing these 
resonances is adopted. We determine the energy of 7?(p) 
as the ground state energy obtained from the correlation 
matrix with ^2 e lpx q{x)Tq{x) interpolators. This energy 
is used for the horizontal lines in the plots. This basis 
renders in all cases just one low-lying state. Within our 
approximation this low-lying state corresponds to a reso¬ 
nance R with momentum p, to a two-particle state mr/nri 
with total momentum p, or to some mixture of R and 
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the two-particle state. We also do not consider nonreso¬ 
nant three-meson levels which could appear above r/ c 7T7r, 
T] c KK, J/ipKK thresholds. Based on the expe¬ 
rience with two-meson operators we do not expect that 
without explicit incorporation of three-meson interpolat¬ 
ing fields these three-meson states appear in the spectra. 

III. LATTICE METHODOLOGY 

These calculations are performed on Nf = 2 dynam¬ 
ical gauge configurations with ~ 266 MeV J28j and 
with other parameters provided in Table HI The mass- 
degenerate u/d quarks are based on a tree-level improved 
Wilson-clover action. The strange quark is present only 
in the valence sector and we assume that the valence 
strange content could uncover hints on the possible ex¬ 
istence of the ccss exotics. The absence of dynamical 
strange quarks prevents ccss intermediate states in the 
cc{uu + dd) and cc sector, in accordance with treating 
these two 1 = 0 sectors separately in our study. With 


a rather small box size of L ~ 2 fm, one expects to have 
large finite size effects. On the other hand this serves as 
a crucial practical advantage by reducing the number of 
two-meson scattering states Mi(n)M 2 (—n) in the energy 
range of our interest. This helps in easier identification 
of the possible resonances that could exist along with the 
regular two-meson energy levels. It also reduces compu¬ 
tational cost as one needs to consider a smaller number 
of distillation eigenvectors and two-meson interpolators 
with respect to a study in larger volume. 


Lattice size 

K P 

N cfgs 

m n [MeV] 

a [fm] 

L [fm] 

16 3 x 32 

0.1283 7.1 

280 

266(3)(3) 

0.1239(13) 

1.98 


TABLE I. Details of the gauge field ensemble used. 

We construct altogether 22 interpolators with J PC = 
1 ++ and total momentum zero for the three cases of our 
interest (T p+ irreducible representation of the discrete 
lattice group Oh is employed): 


Of_ 8 = cMc(O) |(1 + K d ), see Table X of Ref. |H (2) 

= &y 5 u(0) u"fic(0) - cjiu(0) u"f 5 c(0) + K d {u -A d}, 

Oio M = djk C7jc(0) { U7fcM(0) + K d {u -A d} }, 

OT 1 = (^75w(e p ) u'yic(-ep) - c^iu(e p ) M7 5 c(-e p )} + K d {u -A d}, 

e p =±e XtVtZ 

° 12 M = C7574 m( 0) M7i7 4 c(0) - C7i7 4 u(0) U757 4 c( 0) + K d {u -A d}, 

° 13 M = e ijk C7 j 7 4 c( 0) {u7fe7 4 M(0) + K d {u -A d}}, 
qMM _ ^2 £ijl C7jc(e p ) {u~/iu(-e p ) + K d {u -A d}}, 

Cp=ie x ,y,2 

Oif M = {c7 5 c(e p ) uu(—e p ) - C7 5 c(-e p ) uu(e p )} p=l + K d {u -A d}, 

OT 1 = Cij P {cjj 75c(—e p ) u"f 5 u(e p ) - cry^ b c{e p ) v,'y 5 u(-e p )} + K d {u -A d}, 

Or = 07,750(0) uu(0) + K d {u -A d}, 

Oii M = {cc(e p ) M7 5 u(-e p ) - ce(-e p ) u'y*,u{e p )} p =i +K d {u^ d}, 

°19 = [c C l5 unrcu]- 3c + [c Cl rhAc T l,Cu\ 3c + K d {u -a d}, 

°20 = [c Cu P Kira 5 Cu]- 3c + [c C'y i 'y 5 u T ] 3c [c TCu hc + K d{u -A d}, 

°2i = [c C , 75?2 T ] gt: [c T 7 i C'w]6 0 + [c C'7 i 'u T ]gJc T 75C'u] 6c + K d {u -A d}, 

°22 = [c C , u T ]g c [c' r 7 i 75C'w] 6o + [c C'7i75'u T ]g c [c T Cu]e c + K d {u -A d}. 


The indices i, j, k and l define the Euclidean Dirac gamma 
matrices, while the index p indicates the momentum di¬ 
rection. Einstein’s summation convention is implied for 
repeated indices. The unsummed index i in all the op¬ 
erators defines the polarization. The C = ^7274 is the 
charge conjugation matrix. The coefficient K d depends 
on the quark content: K d = 1 is used for cc(uu + dd) and 
K d = 0 for ccss followed by using strange quark propa¬ 
gators instead of the light quark propagators. For 1 = 1 


channel we apply K d = — 1 which gives the flavor content 
cc{uu — dd) and has the same spectrum as ccdu in the 
isospin limit. 


We emphasize the use of four operators 0 4q with 
diquark-antidiquark structure and color antitriplet or 
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N 

cc(uu + dd) 

ccud 

ccss 

nee 

'-'1—8 

c M c 

Does not couple 

c M c 


D(Q)D*(0) 

D(0)D*(0) 

Ds(0)D* s (0) 

n MM 

^10 

j/V/omo) 

j/mp(o ) 

j/mm 

r\MM 

U 11 

D(1)D*(-1) 

D(1)D*(-1) 


'-'12 

D{0)D*{0) 

D{0)D*(0) 

D s (0)D*(0) 

r\MM 

'-'13 

JH>( 0)w(0) 

j/mp(o) 

j/mm 

r^MM 

'-'14 

JMiM-i) 

J/V>(i)p(-1) 

j/mn-i) 

n MM 

'-'15 

r/ c (lM-l) 

??c(l)no(—1) 

Not used 

r-\MM 

'-'16 

Xd(l)»7(-1) 

Xoi(1)tt(-1) 

Not used 

n MM 

(-'17 

Xcl(0)cr(0) 

Xci(0)a o (0) 

Not used 

r-\MM 

(-'18 

Xco(l)f?(—1) 

Xco(1)tt(-1) 

Not used 

o 4q 

^19-20 

[cq]s c [cq\ s„ 

[cm] 3c Msc 

NscNsc 

o 4q 

w 21-22 

[c?] S r . Nfic 

[cm] 6c [cd] 6a 

Msc Mec 


TABLE II. List of interpolators ( J PC = 1 ++ ) and their cor¬ 
respondence with various two-meson scattering channels. 


sextet diquarks 

[crig] s [ c r 2 g] e - = 

Xl 

• E^^C r 2 2/32 4 2 (x2,t/). (3) 

x 2 

Here a = 1, 2, 3 for color triplet and a = 1,..., 6 for sextet, 


while b, c = 1,2,3 for both: 

G a bc = G a bc = e abc (4) 

Gibe = Gl bc = 1 : a = 1,2,3 and a ± b £ c 
Gibe = Gibe = v / 2:a = 4,5,6and«-3 = 6 = c 

while the remaining Q a b c are zero. The operator [Eq. ([3])] 
reduces to c(x)g(x)c(x)g(x) on ensemble averaging, 
where the gauge configurations are not gauge fixed. 

The interpolators are related with the two-meson chan¬ 
nels as listed in Table IE noninteracting levels corre¬ 
sponding to some of these two-meson channels lie above 
our energy of interest, and the corresponding interpola¬ 
tors are not considered. 

The Wick contractions considered in the computation 
of the correlation functions are shown in Figure [T] There 
are two other classes of diagrams, which are not consid¬ 
ered: one in which no valence quark propagates from 
source to sink, and the other class in which only the 
light/strange quarks propagate from source to sink and 
the cc pair annihilates. The effects from these two classes 
of diagrams, with the charm quark not propagating from 
source to the sink, are known to be suppressed due to 
the Okubo-Zweig-Iizuka rule. They correspond to mix¬ 
ing with a number of channels that contain only the u/d 
and s quarks, which represents currently unsolved chal¬ 
lenge in lattice QCD. Note that the annihilation of u/d 
and s quarks as well as mixing with cc is taken into ac¬ 
count, unlike in the simulation [19| aimed at F(4140), for 
example. 



FIG. 1. The Wick contractions considered in our calculations, (a) Connected contraction diagrams, (b) Diagrams, in which 
the light/strange quarks do not propagate from source to sink. The correlation functions in the cc(uu + dd) and ccss cases 
are linear combinations of the diagrams of kind (a) and (b), while the correlation functions between the operators with quark 
content ccud are constructed purely from diagrams of kind (a). 


Using the interpolators listed in Eq. d2J and Table El 
we compute the full coupled correlation functions 

C jk (t) = (niOj^+tplitM =Y J Zrz;e~ E " t . (5) 


For an efficient computation of these correlation matri¬ 
ces, we utilize the “distillation” method for the quark 
sources as proposed in Ref. [29(. In this method the 
quark sources are build from the N v lowest eigenmodes 
of the gauge-covariant Laplacian on a given time slice, t s . 
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We use N v = 64 for computation of correlators involving 
u/d quarks, while for the correlators with hidden strange 
content, we use N v = 48. The correlation functions with 
u/d quarks are computed only for polarization along the 
x-axis and averaged over all t s , while correlation func¬ 
tions involving hidden strange quarks are averaged over 
all polarizations and for all even values of t s . 

The energies E n and overlaps Z = (f2|Oj|n) for 
all eigenstates n are extracted using the well-established 
generalized eigenvalue problem (30] 


C(t)u ( - n \t) = X^(t,t 0 )C(t 0 )u^(t). (6) 

The energies E n are extracted asymptotically from two- 
exponential fits to the eigenvalues 


A (n) (i, t 0 ) (x A n e ~ E + A' n e~ E ^ , E' n > E n . (7) 

We find consistent results for to = 2, 3 and present the 
results for to = 2. The two-exponential fits were typically 
done in the range 3 < t < 14. The eigenvectors u^ 
determine the overlaps 



E„t/2 I Cjk(t)u\, ^(f)| 
\C(t)uW(t)\ 


( 8 ) 


The statistical errors obtained using single-elimination 
jackknife analysis are quoted throughout. 

The complete basis was used in the initial analysis, 
which was later reduced to an optimized basis, separately 
in each of the three cases, based on a systematic opera¬ 
tor pruning. This procedure is aimed at getting better 
signals (in terms of the numbers of states and the qual¬ 
ity of the effective mass plateau and the overlap factors) 
in comparison with the spectrum extracted from the full 
set of operators. After finalizing the optimized set of 
two-meson interpolators, we fixed the cc and [cg]g[cg]q 
operators that give good signals for a maximum number 
of extractable states below 4.2 GeV. The optimized basis 
that we used for the three cases of quark content are 

cc(uu + dd): 01/ 35 , Ogi^ ;14jl5jl7 , 0-^ 2l 
ccud : 0$L% tia , 0\l 21 

cess : Oy 5 , 0]J 21 ■ (9) 

Our principal aim is to find out whether QCD sup¬ 
ports exotic states in addition to the conventional char- 
monia and the two-meson scattering levels, which in¬ 
evitably appear in dynamical QCD. Analytic techniques 
have been proposed for the determination of the scat¬ 
tering matrix for coupled two-hadron scattering chan¬ 
nels based on Liischer-type finite volume formalisms [3l| ■ 
These would in principle allow extraction of the masses 
and decay widths for resonances of interest. A number 
of lattice calculations have already dealt with resonances 
and shallow bound states in the elastic scattering (see 


[32} and [13 for an example of each). The first calcula¬ 
tion of a scattering matrix for two coupled channels also 
promises progress in this direction [34|]. However, such an 
analysis is beyond the scope of current lattice simulations 
for more than two coupled channels and/or three-hadron 
scattering channels, which applies to the case considered. 

Therefore we take a simplified approach, where the ex¬ 
istence of possible exotic states is investigated by ana¬ 
lyzing the number of energy levels, their positions and 
overlaps with the considered lattice operators (fi|0 7 jn). 
The formalism does predict an appearance of a level in 
addition to the (shifted) two-particle levels if there is a 
relatively narrow resonance in one channel. We have, for 
example, found additional levels related to the resonances 
P S3, AT*(892) m, £>5(2400) SI, and the bound state 
DA(2317) S3. Additional levels related to 775(1430) 
SI and A(3872) El have been found in the simulations 
of two coupled channels. Based on this experience, we 
expect an additional energy level if an exotic state is of 
similar origin, i.e. if it corresponds to a pole of the scat¬ 
tering matrix near the physical axis. 

Consider a noninteracting situation. Several two- 
meson operators considered in TablcHTIcontain the vector 
meson V(l) with one unit of momentum. This can reside 
in irreducible representations (irreps) A\ or E 2 of the cor¬ 
responding symmetry group Dic^ [3a [37j. One expects 
two degenerate energy levels for P(1)V(— 1) since there 
are two ways to combine the vector-meson irrep (Aq, E 2 ) 
with the pseudoscalar-meson irrep (A 2 ) to obtain the rest 
frame irrep of interest Tj + (see Table III of [13). The un¬ 
derlying reason is that PV state with J p = 1 + can be 
in s-wave or in d -wave (also in continuum) [38l - l40i |. In 
the limit of small coupling between s— and d— wave, one 
energy level is due solely to the s-wave and the other one 
to d -wave [39|, E3 0 We implement only the s-wave in¬ 
terpolator O p (b l ( -1 ) [Eq. J2J] and therefore expect to 
see only one energy level; this is verified in our observed 
spectra shown in Sect. [V] One would need to employ two 
distinct interpolators in order to find two P(1)V'(—1) en¬ 
ergy levels, but the extraction of such eigenstates has 
not been attempted yet for two-meson systems in QCD 
to our knowledge. Our two-meson operators contain also 
Fi(l)y 2 (—1), where three levels are expected based on 
analogous arguments [13; we expect to find only one 
level related to s-wave interpolators [Eq. <[2j)] , and in¬ 
deed we do not find two other levels related to d -wave 
(for total spins S = 1,2). We emphasize that the omis¬ 
sion of additional interpolator structures and avoidance 
of levels related to d -waves makes the search for possible 
exotics within our approach less cumbersome and results 


1 The PV spectrum resembles (in the noninteracting limit) the 
spectrum in the deuterium channel pn, since 5=1, J p = 1+ 
and l = 0, 2 apply in both cases. Figure 2 of [Tt| indicates that 
one level n(l)p(—1) is related mostly to s-wave and the other to 
d- wave. Liischer’s quantisation condition [Til does not depend 
on the spins of the individual particles, but on their total spin S. 
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more transparent. @ 

Charm quarks being heavy are subject to large dis¬ 
cretization errors. We treat the charm quarks using 
the Fermilab formulation [4lj, according to which we 
tune the charm quark mass by equating the spin aver¬ 
aged kinetic mass of the IS 1 charmonium to its physical 
value. With this formulation, the discretization errors 
are highly suppressed in the energy splitting E n — m s . a ., 
m s .a. = j{m Vc + Smj/i/,), which will be compared with 
the experiments. We utilized this method in our ear¬ 
lier calculations on this ensemble and found good agree¬ 
ment with the experiments for conventional charmonium 
in Ref. [25J as well as, for masses and widths of charmed 
mesons in Refs. !25 ; , 33i. H3|. 


IV. FIERZ RELATIONS 


The diquark-antidiquark operators [cq\s a \cq\^ c and 
[cg]g c [c(7]6 c can be expressed as linear combinations of 
color singlet currents (cc)i c (qq)i c and (cq)i c {qc)i c [H|, 
|43| . These relations are obtained for local currents via 
Fierz rearrangement [44| and are presented in the Ap¬ 
pendix. Note that our quarks are smeared and each 
meson in O mm has definite momentum, but the Fierz 
relation suggests that O iq and O mm are still linearly 
dependent. 

The Fierz rearrangement is the key idea behind Cole¬ 
man’s argument [45j] that in the large N c limit application 
of Fermion quadrilinears to the vacuum creates meson 
pairs and nothing else. In the physical world with N c = 3, 
it is argued that tetraquarks could exist at subleading or¬ 
ders [4|| of large N c QCD. However, in the presence of 
the leading order two-meson terms, one should take cau¬ 
tion in interpreting the nature of the levels purely based 
on their overlap factors onto various four-quark interpo¬ 
lators. 

Let us consider a comparative study between the lat¬ 
tice correlators and the Fierz expansion of O iq operators. 
From Eq. (|A5j) . we see that the first and second terms 
in the Fierz expansion represent DD *, while the seventh 
term is similar to the Of = Xci a. Hence we expect 
significant correlations between these operators. This is 
indeed verified in Figure (2[ showing the time averaged 
normalized ensemble averaged correlation matrix 


!J 9 hi y/CaiWnit) 


( 10 ) 


With this normalization all the diagonal entries are 
forced to unity and all the off-diagonal entries to be less 


2 If one would find an extra state near V(1)P(— 1) or Vi (1) V 2 (—1), 
one would indeed have to identify whether this extra state arises 
due to the presence of the d- wave or is related to exotics. We 
do not address this question since we do not find such an extra 
state. 



FIG. 2. Time averaged normalized correlation matrix C [Eq. 
pin for the operator basis Oi _22 with quark content cc(mm + 
dd) and cc. The axis ticks correspond to the order of operators 
used in Eq. ©• 


than unity. The [cg]g[cg]e = O 19-22 have large correla¬ 
tions onto the DD* = Ogf^ 12 and Xci 7 = Off M . The 
strong correlations between 0 4q and O cc operators can 
also be explained by the Xcic component in 0 4q : where 
a couples to the vacuum. 


V. RESULTS 


The discrete spectra in Figs. [3] and [4] are the main 
results from our lattice calculation. They show the ener¬ 
gies 


E n = El at - m‘ 


~lat 


m 


exp 


• 3777-jy^) 


(ii) 

of the states with J PC = 1 ++ and three quark contents. 
The horizontal lines represent various two-meson nonin¬ 
teracting energies. 

The states that have dominant overlap with two-meson 
scattering operators are represented by circles and the 
color coding identifies the respective scattering channels 
based on the following criteria: 


• The levels appear close to the expected two-meson 
noninteracting energies. 

• They have dominant overlaps (fl\Of IlAl2 \n) with 

corresponding This is also verified based 

on the ratios /max m {ZJ l ) 1 which are indepen¬ 
dent of normalization of operators and are shown 
in Figure [6[ 


• If the corresponding two-meson interpolators are 
excluded from the basis, this eigenstate disappears 
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(a) 1=0 : cc ( uu + dd ) & cc 


JW(1)<B(-1) 


D(1)D* (-1) 

Xci(°) o(0) 

JW(0J <o(0) 




Zcidl^M) 

Wb^-I) 

D(1)D* (-1) 
JAt-(l) p(-1) 

J/'i'tO) p(0) 
D(0) D*(0) 


FIG. 3. The spectra of states with J PC = 1 ++ for the cases with u/d valence quarks. The energies E n = E ^ ai — rris.a. + mt x £. 
[Eq. (Illll l are shown. The horizontal lines show energies of noninteracting two-particle states JTJ and experimental thresholds, 
indicating uncertainty related to a width. In each subplot, the middle block shows the discrete spectrum determined from our 
lattice simulation from the optimized basis [Eq. (|9))]. The right-hand block shows the spectrum we obtained from the optimized 
basis of operators with the [ cq]g[cq]g operators excluded. The gray marks, on the right-hand side of each pane, indicate the 
lowest three-meson threshold m^ c + 2 m n , while the actual lowest r/ c nn level on the lattice appears higher due to l = 1, which 
requires relative momenta. The left-hand block shows the physical thresholds and possible experimental candidates (a) Xci, 
X(3872) and A'(3940), (b) .£+(4050) and .£+(4250). The violet error bars for experimental candidates show the uncertainties 
in the energy and the black error bars show its width. 
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FIG. 4. The spectrum of states with J pc = 1 ++ and hidden 
strange quarks. The possible experimental candidates shown 
are Xci, A(3872), T(4140) and V(4274). The gray marks, on 
the right-hand side of each pane, indicate the lowest three- 
meson threshold m Vc + 2m,K- However, the actual lowest 
rf c KK level on the lattice appears higher due to l = 1, which 
requires relative momenta. For further details see Figure [3] 


or becomes too noisy to be identified. This is de¬ 
termined by comparing the pattern of the effective 
masses and overlaps between the original basis and 
the basis after operator exclusion. 

The remaining states, that are not attributed to the 
two-meson scattering channels, are represented by red 


FIG. 5. The spectrum of states (Eq. (1111) 1 with J PC = 1 ++ 
and quark content cc(uu + dd) & cc. (i) Optimized basis 
(without Oi 7 M ), (ii) optimized basis without cc operators 
(and without Oiy M ) and (iii) basis with only cc operators. 
Note that candidate for A'(3872) disappears when remov¬ 
ing cc operators although diquark-antidiquark operators are 
present in the basis, while it is not clear to infer on the dom¬ 
inant nature of this state just from the third panel. The 
Of^ M = Xci(0)<r(0) is excluded from the basis to achieve bet¬ 
ter signals and clear comparison. 


squares. 

Figures [3] and 0] also compare the spectra between the 
two bases of operators, one with optimized operator set 
and another with the optimized set excluding [cq]g[cq]g. 
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FIG. 6. The overlap factors Z^ n) [Eq. ([8])] shown in units of the maximal | Z™ \ for a given operator j across all 

the eigenstates m . These ratios are independent of the normalization of the interpolators Oj. The horizontal axis corresponds 
to the complete basis of interpolators [Eq. (|2)l]. where the optimized subsets [Eq. ([9])] were employed. The levels are ordered 
from lowest to highest E n as in the middle pane of the spectrum in Figs. [3]and[I| The values are averages of the ratios over 
4 < / < 13 with error bars due to jackknife sampling. 
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In all three cases we see an almost negligible effect on the 
low lying states, while we do observe an improvement in 
the signals for higher lying states in the basis without 
[i cq]g[cq]g. The same conclusion applies for overlaps. 

The employed irreducible representation Tj ++ con¬ 
tains the states J PG = 1 ++ of interest, as well as 
J PC = 3 ++ states due to the broken rotational sym¬ 
metry. Upon inclusion of the interpolator 0§ c to the 
basis [Eq. @] the spectra for both 7 = 0 channels re¬ 
main essentially unchanged except for an additional level 
at E ~ 4.1 — 4.2 GeV [Eq. (fill) ]. This is where the earlier 
simulation on the same ensemble [25[ and the simula¬ 
tion [47] have identified the only 3 ++ state in the energy 
region of our interest. In the following subsections, we 
present the spectra of J PC = 1 ++ states in three flavor 
channels for the basis (Eq. @), where 0§ c is excluded. 


A. 7 = 0 channel with flavor cc(uu + dd) and cc 

This is the channel where the experimental A (3872) 
resides. We will argue that the energy levels affected 
by this state are n = 2 (red squares) and n = 6 (blue 
circle) from Figure [3][a). The lowest state is the con¬ 
ventional Xci(l-P)- The overlaps of the three low-lying 
levels represented by circles show dominant J/^(0)w(0), 
77c(l)er(—1) and Xci(0)cr(0) Fock components. The high¬ 
est two states in Figure[3][a) have significant overlap with 
the J/ip(l)uj(—l) and D 0 (1)Dq(— 1) operators. 

Now we focus on the eigenstates that are related to 
A (3872). The cc interpolators alone give an eigenstate 
close to DD* threshold (right pane of Figure O, but 
one cannot establish whether this eigenstate is related 
to A(3872) or to nearby two-meson states in this case. 
Therefore we turn to the spectrum of the full optimized 
basis [midpane in Figure [3](a)], where levels n = 2 (red 
squares) and n = 6 (blue circles) are found to have dom¬ 
inant overlap with the cc and DD* operators. Exclud¬ 
ing either of these operators results in disappearance of 
one level and a shift in the other level towards the DD* 
threshold. We emphasize that one of the two levels re¬ 
mains absent when DD* and 0 Aq are used and O cc is not, 
as is evident from the first and second panel from the left 
of Figure [5] This indicates that the cc Fock component is 
crucial for A(3872), while the [cq\g[cq]g structure alone 
does not render it. This also implies a combined domi¬ 
nance of cc and DD* operators in determining the posi¬ 
tion of these two levels, while their resulting energies are 
not significantly affected whether 0 Aq is used in addition 
or not. 

We determine the DD* scattering phase shift from lev¬ 
els n = 2, 6 via Liischer’s relation [3JJ assuming elastic 
scattering. The phase shift is interpolated near threshold 
using the effective-range approximation. The eigenstate 
n = 6 (blue circle) is interpreted as the 71(0)7)* (0) scat¬ 
tering state, which is significantly shifted up due to a 
large negative scattering length j48[. The resulting scat¬ 
tering matrix T oc l/(cot<5(p) — i) has a pole just below 


A(3872) 

mx — rris.a. 

mx — mu, - m D * 

Lat. 

816(15) 

-8(15) 

Lat. - 0 4q 

815(8) 

-9(8) 

LQCD [17] 

815(7) 

-11(7) 

LQCD [18] 

- 

-13(6) 

Exp. 

803(1) 

-0.11(21) 


TABLE III. Mass of A'(3872) with respect to m 3 . a . and the 
DqDq threshold. Our estimates are from the correlated fits 
to the corresponding eigenvalues using single exponential fit 
form with and without diquark-antidiquark operators. Re¬ 
sults from previous lattice QCD simulations [in . |T§] and ex¬ 
periment are also presented. 


the threshold where cot 8{pb) = i is satisfied. We neglect 
possible effects of the left-hand cut in the partial wave 
amplitude. The results confirm a shallow bound state 
just below the DD* threshold and the binding momen¬ 
tum pb renders the mass of the bound state, interpreted 
as experimentally observed A (3872). The resulting mass 
of A(3872) and its binding energy are provided in Table 
uni and in Figure [3 which indicate that it is insensitive 
to inclusion of diquark-antidiquark interpolators within 
errors. The mass of A(3872) was extracted along these 
lines for the first time in Ref. [l7j, where this channel 
was studied in a smaller energy range on the same en¬ 
semble without diquark-antidiquark interpolators. The 
error on the binding energy in the present paper is larger 
due to the larger interpolator basis. These results are 
in agreement with a possible interpretation of X(3872), 
where its properties are due to the accidental alignment 
of a cc state with the D°D*° threshold [49|, j50j, but we 
cannot rule out other options. 


850 


790 

770 

10 

0 


-30 


m X(3872) m s.c 


♦ 

+ ♦ 

m 


m X(3872)~ m D 

-m~ D * 

- ♦ 

! ♦ 

■t it ; 

Exp. 

Lat. Lat .-0 4q 

[17] [18] 


FIG. 7. Mass of A'(3872) with respect to m s . a . from the 
present simulation, previous lattice studies [Til Fl8 ] and ex¬ 
periment [fl]. 

With regard to the other experimentally observed 
charmonia-like states [e.g. A(3940)], which could appear 
in this channel, we do not find any candidate in addition 
to the expected two-meson scattering levels. We also do 
not find candidates for other cc states with J PC = 1 ++ 
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[e.g. xi c(nP)\ in the region between the DD* threshold 
and 4.2 GeV. 


B. 1=1 channel with flavor ccdu 

A careful analysis of this isospin channel is crucial due 
to the large branching ratio for the decay V(3872) —> 
J/ipp and current experimental interests in search of a 
charged A! (3872). With no disconnected diagrams al¬ 
lowed in the light quark propagation, the correlation ma¬ 
trix is constructed purely of four-quark interpolators and 
connected Wick contractions in Figure (Ha). 

The spectrum of eigenstates is shown in Figure [3jb) , 
where all energies are close to noninteracting energy lev¬ 
els. All the eigenstates have a dominant overlap with the 
two-meson interpolators. The spectrum shows very little 
influence on the inclusion of [cgjgfcgjg, which is evident 
from Figure02b). Given that all the levels below 4.2 GeV 
can be attributed to the expected two-meson scattering 
states, we conclude that our lattice simulation gives no 
evidence for £ c (4050)+ and £ c (4250)+. 

Our results also do not support charged or neutral 
V(3872) with 1 = 1 . There is no experimental indi¬ 
cation for charged X , while the neutral X does have a 
large decay rate to I = 1 final state J/V'p 0 - One popu¬ 
lar phenomenological explanation for this decay is that 
V(3872) has 1 = 0 and the isospin is broken in the decay 
mechanism (due to the D + D*~ vs D°D*° mass differ¬ 
ence) [50, |5l(. According to another explanation, X is 
a linear combination of I = 0 and 1 = 1 components, 
where the 1 = 1 component vanishes in the isospin limit 
[ 52 }. Our simulation is performed in the isospin limit 
m u = rrid : so it is perhaps not surprising that X with 
1 = 1 is not observed. Future simulations with nondegen¬ 
erate u/d quarks would be very welcome for this channel. 

As pointed out in Section m p in J/V^p is treated as 
stable, although p{ 1) is kinematically close to the de¬ 
cay channel 7r(l)7r(0). In the absence of a simulation of 
a three-meson system, it is disputable what ‘noninter¬ 
acting’ energy should be taken for the p( 1). An esti¬ 
mate from the diagonal correlator p(l) leads to ‘nonin¬ 
teracting’ energy roughly 65 MeV below the eigenstate 
energy, which is identified to have a dominant overlap 
with the J/ip(l)p{— 1) interpolator. However, taking the 
resonance position [261] brings the ‘noninteracting’ level 
in agreement with the measured eigenenergy. 


C. 1 = 0 channel with flavor ccss and cc 

Our goal in simulating this channel is to search for a 
possible presence of the Y(4140) resonance, which was 
found in J/4’ ( t ) scattering in several experiments 0- 
Our lattice simulation of J/ip(j> scattering takes into ac¬ 
count the annihilation of the valence strange quarks and 
thereby the mixing with cc flavor content. 


With no strange quark effects in the sea, the study 
of this channel is based on the following assumptions. 
We construct a basis with only cc and four-quark oper¬ 
ators (O mm , O iq ) with valence hidden strange content 
for this analysis. We assume that these interpolators 
have negligible coupling to two-meson states with flavor 
content c.c(uu + dd). In other words, we assume that 
two-meson states like DD* and J/ipco will not appear 
in the spectrum based on the chosen interpolators. The 
resulting spectrum in this channel confirms this assump¬ 
tion. We point out that Y(4140) has been experimen¬ 
tally observed only in the final state with valence 

strange content, but it has not been observed in DD* and 
J/ijjtjU final states. Although this ensemble does not have 
strange quarks in the sea, we assume that the valence 
strange content could uncover hints on the existence of 
the charm-strange exotics, if they exists. 

Spectra in this channel are shown in Figure |U We 
identify the lowest two states, represented by squares, 
to be Xci(l-P) and the level related to X(3872). The 
remaining four states are identified with the expected 
D S D* and scattering levels. Thus in the energy 

region below 4.2 GeV, we find no levels that could be 
related to Y(4140) or any other exotic structure. Note 
that the existence of Y(4140) is not yet finally settled 
from experiment, and its quantum numbers, except for 
C = +1, are unknown. Therefore it is possible that its 
absence in our simulation is related to the fact that we 
explored the channel J p = 1 + only. 


D. Discussion 

The only exotic charmonium-like state found in our 
simulation is a X(3872) candidate with J PC = 1++ and 
1 = 0. It is found as a bound state slightly below DD* 
threshold and has a mass close to the experimental mass 
of V(3872). We point out that this mass corresponds 
to our m,r ~ 266 MeV and was obtained from a rather 
small lattice volume, while chiral and continuum extrap¬ 
olations have not been performed. Precision determina¬ 
tion of its mass with respect to DD* threshold will be 
a challenging task for future lattice simulation on larger 
volumes, which also should account for its coupling with 
multiple open scattering channels involving two or more 
hadrons. Recent analytic studies consider the quark mass 
dependence, the volume dependence and the effect from 
the isospin breaking relevant for future lattice studies of 
X(3872) HI. 

Candidates for no other “exotic” charmonium-like 
states [except for X(3872)] are found in our exploration 
of the three J PC = 1 ++ channels. We list several possi¬ 
ble reasons for the absence of the energy levels related to 
other possible exotic states in our simulation: 

• The existence of Y(4140), Z+{ 4050), £+(4250) or 
any other exotic state in these channels, is not yet 
settled experimentally. Even if they exist, only C = 
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+1 is established experimentally, while their J p is 
unknown. This could explain their absence in our 
simulation, which probes only J p = 1 + . 

• Based on the experience, discussed in Section m 
we expect an additional energy level if the exotic 
state is a resonance associated to a pole near the 
real axis in the unphysical Riemann sheet. The 
absence of an additional energy level could also in¬ 
dicate a different origin of the experimental peak, 
e.g., a coupled-channel threshold effect. Further 
analytical work and lattice simulations are needed 
to settle the question whether an additional energy 
level is expected in this case. 

• Finally, we cannot exclude the possibility that some 
exotic candidates could be absent due to the rela¬ 
tively heavy pion mass m n — 266 MeV, isospin 
limit m u = neglect of the charm annihilation 
contributions, or the absence of the strange dynam¬ 
ical quarks in our simulation. 

VI. CONCLUSIONS 

We present the spectra from a lattice QCD simulation 
of J PG = 1 ++ channels with three different quark con¬ 
tents: ccdu , cc{uu+dd) and ccss , where the later two can 
mix with cc. The pion mass in this study with u/d dy¬ 
namical quarks is m^ ~ 266 MeV. Using a large number of 
interpolating fields^ [cq] sjcg]^, [cq]e c [cq\ 6c , (cq)i c (qc) lc , 
(cc)i c (qq)i c and (cc)i c , we extract the spectra up to 
4.2 GeV. We find evidence for and X(3872), while all 
the remaining eigenstates are related to the expected two- 
meson scattering channels, which inevitably appear in 
the dynamical QCD. The cc Fock component in A(3872) 
appears to be more important than the \cq\g{cq\g, since 
we find a candidate for X(3872) only when cc interpolat¬ 
ing fields are used. The DD* interpolators show a more 
prominent effect on the position of A(3872) than the 
[cq}g[cq]g. Candidates for charged or neutral A(3872) 
with 7 = 1 are not found in our simulation with m u =rrid , 
and future simulations with broken isospin would be wel¬ 
come for this channel. We also do not find a candidate for 
Y (4140) or any other exotic charmonium-like structure. 
Our search for the exotic states assumes an appearance of 
an additional energy eigenstate on the lattice, which is a 
typical manifestation for conventional hadrons. Further 
analytic work is needed to establish whether this working 
assumption applies also for several coupled channels and 
all exotic structures of interest. 
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Appendix A: Fierz transformation of 
diquark-antidiquark operators 

In this appendix we express the local diquark- 
antidiquark interpolator as 

O iq {x) M i( x ) M ^ x ) ( A1 ) 

using Fierz transformations [44j] . where 

M(x) = q a Tq' a (x) (A2) 

are local color-singlet currents. The momentum pro¬ 
jected interpolator 0 4q (p) is then given by 

0 4q (p) = E F i E e ipx Mi(x)M*(x) 

i x 

= E $ E e<P * E e ~ iqX M\{q ) E e~ ikx Mi(k) 

i x q k 

i q 

Thus the projection to total momentum zero 0 4q (p = 0) 
can be rewritten as sum over two-meson operators with 
back-to-back momenta. 

Fierz transformation is an operation of rearranging the 
Fermion fields in a Fermion quadrilinear. Expressing our 
local diquark-antidiquark interpolator with explicit color 
(lower) indices and Dirac (upper) indices, we have 

[C P q]g[ C N g ] S | ( 3e) = QabcGade C$t q* e 

= (SbdSce T SbeScd) P aP N q& Clt C V d q S e 

= P a 0 N q5 {-(c^k(«?V) lc T (cV)i c (?V) lc } 

= -(c“r?V) le (g*G?Y) lo =F (^Tf s ^) lo (fH^) lc 
= -(c Tj c)i c (q Gi q)i c =F (c Tj q)i c {q 77/ c)i c (A3) 

where we have accounted for GabcGade = Sbd$ce T 
Shewed in the second line and a minus sign for Fermion 
exchange in the third. Each term on the right-hand side 
of the fourth line is expressed as a sum over the index 
7 = 1,..., 16, where Tj are the elements of Clifford alge¬ 
bra {r} and (G/, 77/) the unknown coefficient matrices. 
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These coefficient matrices can be determined using the 
orthogonality relation TrpjTj] = ASjj 

Gi = ^(N t TiP) t and id/ = i(lVr z P) T (A4) 

I 


The diquark-antidiquark fields can therefore be ex¬ 
pressed as a linear combination of products of two color 
singlet currents with various Dirac structures. Local 
analogs of our diquark-antidiquark interpolating fields 
can be expressed as 


= [c C 75 u\g[c 7iC u\g + [c C 7i u]g[c 7 5 C u)g + K d {u -7- d} (A5) 

\2l) 

(- 1 )' 

= T 2—{ ( C 75 7 * c ) - ( c liu){u 75 c) 

+ (c 7^75 m)(m 7i7„ c)|j^ - (c 7i7„ u)(zt 7^75 c)|i^^} 

(- 1 )* 

H-— { (c c)(u 7,75 u) + (c 7,75 c)(u u) 

- (c 7"c)(u 7,7^75 u)|i^ - (c cr a/3 c)(m cr a/3 7*75 M)| i7 q a</3) } 

+ K d {u —> d} 


and 



= [c C u]g[c 7 ^ 75(7 u]g + [c C 7 i 7 5 u]q[c C u]g + K d {u ->■ d} 

(- 1 )* 

= T— 2 —{ ~ ( c M )( M 7*75 c) + (c 7,75 m)(m c) 

- (C 7 I/ m)(m 7i7^7 5 C)|j^ + (c (J Q/3 u)(u Cr a/ 37 i 75 c)|^ (a</3) } 
(-l ) 1 2 3 4 5 6 7 8 

-— { (c c)(m 7,75 m) - (c 7,75 c)(u m) 

+ (c 7 I/ c)(m 7,7^75 u)|i^ - (c cr Q/3 c)(m 0-0,^37^75 M)|i^(a<^)} 

+ K d {u —1 d}. 


(A6) 


Various terms resemble two-meson operators O mm 


[Eq. |2], where (glV) denote color singlet currents. 
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